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ABSTRACT 

We study the spiral arm influence on the solar neighbourhood stellar kinematics. As the 
nature of the Milky Way (MW) spiral arms is not completely determined, we study two mod- 
els: the Tight- Winding Approximation (TWA) model, which represents a local approximation, 
and a model with self-consistent material arms named PERL AS. This is a mass distribution 
with more abrupt gravitational forces. We perform test particle simulations after tuning the 
two models to the observational range for the MW spiral arm properties. We explore the ef- 
fects of the arm properties and find that a significant region of the allowed parameter space 
favours the appearance of kinematic groups. The velocity distribution is mostly sensitive to 
the relative spiral arm phase and pattern speed. In all cases the arms induce strong kinematic 
imprints for pattern speeds around 17kms~^kpc~^ (close to the 4:1 inner resonance) but 
no substructure is induced close to corotation. The groups change significantly if one moves 
only ~ 0.6 kpc in galactocentric radius, but ~ 2 kpc in azimuth. The appearance time of each 
group is different, ranging from to more than 1 Gyr. Recent spiral arms can produce strong 
kinematic structures. The stellar response to the two potential models is significantly different 
near the Sun, both in density and kinematics. The PERL AS model triggers more substructure 
for a larger range of pattern speed values. The kinematic groups can be used to reduce the 
current uncertainty about the MW spiral structure and to test whether this follows the TWA. 
However, groups such as the observed ones in the solar vicinity can be reproduced by different 
parameter combinations. Data from velocity distributions at larger distances are needed for a 
definitive constraint. 

Key words: Galaxy: kinematics and dynamics - solar neighbourhood - Galaxy: evolution - 
Galaxy: structure - Galaxy: disc - galaxies: spiral 



1 INTRODUCTION 

The spiral arms of our Galaxy have been typically characterised 
by radio observations of the 21 -cm line of neutral hydrogen, 
giant HII regions, CO emission, an d optical data of young O 



and B stars (e.g. Port, Kerr & W esterhoutlll958|; ISimonsonI [19701 : 
lOeorgelin & GeorgelinI Il976h . [brimmel & Spergdl fcOOlh used 
FIR and NIR emission from the COBE/DIRBE to fit a model 
for th e spiral arms in the st ellar comp onent of the Galaxy. Re- 
cently, ^njaminetaD JIool) and lChurchwell all (l2009l), based on 
infrare d data from the Spitzer/GLIMPSE survey, and iReid et all 
(l2009h . using trigonometric parallaxes and proper motions of 
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masers in high-mass star-forming regions, reported new results 
about the spiral structure of the Milky Way (MW). Despite the ef- 
fort, there are still several caveats in the spiral arm properties in our 
Galaxy, such as pattern speed, strength, orientation and even the 
number of arms or their stellar or gaseous structure. Furthermore, 
the nature of the spiral arms themselves, i.e. their origin or their 
lifetime, are nowadays a matter of debate (Sellwood 2010b). 

Apart from direct methods to detect spiral arm over-densities, 
an alternative method is based on the analysis of kinematic groups 
that are induced by the spiral arms in the local velocity distribu- 
tion. Some of the moving groups in the solar neighbourhood were 
originally thought to be remnants of disrupted disc stellar clusters. 
However, the age or metallicity distribution of their stars are in 
contradiction with this hypothesis (iBensbv et al.l2007l : lAntoia et al.l 
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l2008h . lKalnai3 ( ll99lh suggested that a moving group could be a 
group of stars that crosses the solar neighbourhood following a cer- 
tain type of orbit induced by the Galactic bar gravitational potential. 
The demonstration that the bar resonances can produce a kinematic 
group, similar to the observed Hercules group dPehnenlbOOOl; IFuxI 
l200lk has lead to the exploration of si milar effects on the velocity 
distribution due to the spiral arms (see lAntoja et"aDl2010l for a re- 
view). This type of kinematic structures could depend strongly on 
some characteristics of the bar and spiral arms and, therefore, they 
result useful for our understanding of the MW large-scale structure 
and dynamics in its present and past form. 

Few studies, however, have focused on the spiral arm ef- 
fects on the d isc velocity distribution. By integration of test 
particle orbits, IPe Simone. Wu & Tremainel (Eo04) showed that 
few, but intense stochastic spiral density waves produce struc- 
tures which are arranged in branches, resembling thos e from 
observations (Skuljan, Heamshaw & Cottrell 1999; Antoi a et aP 
boOSh . lbuillen & MinchevI (120051) developed a method to quantify, 
through orbital integration, the velocity distribution that would re- 
sult from the existence of spiral-induced families of periodic orbits 
in the solar neighbourhood. They found, for models with the Sun 
in the outer limits of the 4:1 resonance and certain spiral arm ori- 
entation, two periodic orbits in the kinematic positions of Coma 
Berenices and Hya des-Pleiades. The test particle simulations by 
IChakrabartvl(l2007h showed that the combined effect of a bar and 4 
weak spiral arms, is necessary to reproduce the main local moving 
groups (Hercules, Hyades, Pleiades, Coma Berenices and Sirius). 
In the spiral-only models, more structures but l ess bold th a n those 
with the bar simulations were obtained. Finally, ' SellwoodI (l2010ah 
showed that the angle-action variables of the stars in the Hyades 
stream could be consistent with the effects of a recent inner Lind- 
blad resonance of a multi-arm (number of arms m > 2) and tran- 
sient pattern. 

All this work has shown that the spiral arms are able to induce 
kinematic groups in the local velocity distribution. However, there 
is no study that explores the influence of each spiral arm property 
using potential models designed according to the latest observa- 
tional evidences for the MW spiral structure. On top of that, pre- 
vious studies have modelled the spiral arm potential with a cosine 
function, following the Ti ght- Winding Approximation (TWA, e.g. 
iBinnev & Tremainel l2008h . However, it is not clear to us whether 
the MW spiral arms satisfy the conditions for self-consistency of 
this model, in particular, regarding their tightness and weakness. 
Moreover, the TWA implies that the spiral arms are a steady global 
mode of the disc but, as mentioned, it is still a matter of discussion 
whether the spiral arms are long-lived or transient structures. 

Motivated by the recent picture of the MW spiral arms and 
the possibility that the TWA is not suitable for our Galaxy, we 
study the spiral arm effects on the solar vicinity kinematics with 
the TWA but also a d ifferent model, namely the PERL AS model 
jPichardo et al.ll200"3l) . This is a 3D mass distribution from which 
more abrupt gravitational potential and the forces are derived. As 
proved in Franco et al. (2002), this difference has far-reaching con- 
sequences on the gaseous dynamical behaviour. Differences in the 
stellar response are expected too. In Antoia et al.i (i2009i) we pre- 
sented a preliminary study of the kinematic effects of this model. 
Here we aim to determine the conditions that favour the appear- 
ance of kinematic groups such as the ones that we observe near the 
Sun. We explore the influence of the spiral properties and discuss 
to which extent we can use the kinematic imprints to constrain the 
spiral arm properties and nature. We also compare the local kine- 
matic imprints of the PERLAS model with the ones of the TWA 



to see whether these imprints are useful to test the spiral arm grav- 
itational potential modelling. To do this, we first tune these two 
models to the latest observational determinations for the properties 
of the MW spiral arms. For some properties only a range of evi- 
dence is available in the literature. Then we perform test particle 
simulations with both models considering these ranges, different 
initial conditions and integration times. 

Section|2]deals with the observed properties of the MW spiral 
structure. In Section |3] we elaborate on the gravitational potential 
modelling of the spiral arms. Section |4]presents the two spiral arm 
models and we fit them to the MW spiral arms. In Section |5] we 
contrast the force fields of the two models. Section |6]describes the 
test particle simulations and initial conditions. Afterwards, in Sec- 
tion |7] we compare the response in density and kinematics to the 
two models. We explore the influence of spiral arm characteristics, 
such as pattern speed, orientation and integration time, on the lo- 
cal velocity distribution in Sections [8l|9] and [TO] In Section [TT] we 
discuss whether it is currently possible to constrain the spiral arm 
properties by using the local observed kinematic groups. Finally, 
we summarise our main results and their implications. 



2 THE SPIRAL STRUCTURE OF THE MW 

Most of the properties of the MW spiral structure remain rather 
undetermined. Here we establish a range of plausibility for each 
property in order to tune our spiral arm models to the MW spiral 
arms. Table [T] shows these adopted ranges. 

Locus. The geometry of the MW spiral arms is still a matter 
of intense debate. Different estimates of the pitch angle and even 
the nu mber of arms can be found from different tracers (see IVallej 
l2008h . Maps of OB -associations and Hll -regions, CO emission, 
mase r s in high-mass star-forming regions (iGeorgelin & GeorgelinI 
"1976", "Tavlor & Cordesl ll993L IValleell2008l and references therein, 
Reid et al. 2009) show a 4- armed pattern, usually referred as 
Sagittarius-Carina, Scutum-Centaurus (or Scutum-Crux), Perseus 
and the outer (or C ygnus) arms. However, according to COBE K- 
band observations (IDrimm el & Spergel 2001) and to the infrared 
Spitzer/GLIMPSE survey (iChurchwell al.. 2009). only two of these 
four arms are major MW arms. These are traced by the stellar 
(young and mainly old) population as enhancements in the spi- 
ral tangencies. External galaxies often show different morpholo- 
gies in blue and near-infrared colours (Grosb0l, Patsis & Pom- 
pei 2004). Although it might not be a general property of galax- 
ies ( Eskridge et al. 2002), many galaxies classified as flocculent or 
multi-armed syste ms in blue, display a 2-armed grand design spi- 
ral in the K band (iBlock etalJll994 Kendall, Kennicutt & Clarke 
201 1). Also hydrodynamic models have shown that it is possible to 
form arms of compressed gas, without increasing the stellar surfac e 
density, as a response to a 2-armed pattern ( Martos et al ]|200i). 
Here we consider only these two stellar major arms (m = 2) as 
they trace the underlying mass distribution of our Galaxy. 

The exact location of the stellar spiral arms is still unclear and 
part of the current models in the outer regions of the MW are an 
extrapo lation of observations in th e inner region s. iBeniamin et al.l 
(2005) JChurchwell al.l (l2009h and I Vailed (l2008h consider that the 
stellar arm counterpart of Scutum-Centa urus is the Perseus arm. 
But according to the model by Drimme Tcfe Spergei| (|200l|) , the 
counterpart goes far beyond the Perseus arm in the anti-centre di- 
rection. Because of these discrepancies, we wil l adop t two different 
locus: the one reported bv lDrimmel & SpergellfcOOlh (hereafter lo- 
cus 1 , solid black curve in Fig.[T]), and the fitting bv I Vailed (l2008h 
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Figure 1. Adopted locus for our spiral arm models: locus 1 (solid line) and 
locus 2 (dashed line). The star at X = and Y = 8.5 kpc indicates the 
assumed solar position. Open circles indicate several regions near the solar 
neighbourhood. 



Table 1. Assumed values or ranges for the properties of the MW spiral arms. 



Property 




Value or range 


Number of arms 


m 


2 


Scale length 


i?s(kpc) 


2.5 


Locus beginning 


Rsp ( kpc) 


2.6/3.6 


Pitch angle 


in 


15.5/12.8 


Relative spiral phase 




88/60 


Pattern speed 


Qsp (kms^-'^kpc"-'^) 


15-30 


Density contrast 


A2 


0.14-0.23 


Density contrast 


K 


1.32-1.6 



for the Scutum and Perseus arms (hereafter locus 2, dashed black 
curve in Fig.[T]). 

We use the galactocentric cylindrical coordinates {R, cj)) with 
the azimuth > in the direction of rotation and origin as indi- 
cated in Fig. [T] The locus of one arm is obtained through solving 
2(0 — 00) -\-g{R) = 0, with the condition ^ 00. The second arm 
is obtained by symmetry. The constant 0o fixes the arm orientation 
and it corresponds to the azimuth of the line that joins the two start- 
ing points of the spiral locus. It is —20° for locus 1 and —70° for 
locus 2. The function q(R) defines the spiral shape. We adopt the 
one from [Roberts. Huntlev & van Albadal(ll979l) : 




where Rsp is the radius of the beginning of the spiral shape locus. 
The parameter N measures how sharply the change from a bar-like 
to spiral-like occurs in the inner regions. Here the limit of ^ oo 
is taken, approximated by = 100, which produces spiral arms 
that begin forming an angle of ^ 90° with the line that joins the two 
starting points of the locus. We have checked that our simulations 
(which are studied in the outer disc) do not depend on the exact 
shape in these inner parts. 

The pitc h angle is i = l^-^. .^ ^-^ "z locus 1 and 2, as 
estimated by Drimmel & SpergellfeOOlh and lValleel(l2008h . respec- 
tively. The relative spiral phase (I)sp{Rq) is the azimuth between 
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Table 2. Density contrast of the MW spiral arms from several studies. 



Author 


K 


A2 


Drimmel & Spergel (2001) 


1.32 


0.14 


Beniamin et al. (2005) 


1.30 


0.13 


Grosb0l et al. (2004) 


1.2-1.6 


0.1-0.23 



the Sun's position and the peak of the spiral at the same radius 
(curved long arrow in Fig.[D- It is 88° and 60° for locus 1 and 2, re- 
spectively. The Sun's is initially assumed to be at = 8.5 kpc and 
= 0° (star in Fig. (TJ. But we will examine also nearby regions 
(open circles in Fig.[T]) for which the relative spiral phase varies ap- 
proximately from —45° to 45° for locus 1, and from —85° to 10° 
for locus 2. With this we are more flexible in the spiral arm location 
which, as mentioned, is rather undetermined. 

Pattern speed. Estimations for the spiral structure pattern 
speed Qsp come from open clusters birthplace analysis, kinematics 
of young stars, and comparisons of t he observed ^ ^CO (l,v) dia- 
gram with models for the gas flow (see 'Gerhard"20 1 for a re view) . 
For the two- armed K-band model by Drimmel & Spergdl (l200lh 
used in the present study, a pattern speed of 20 km s~^kpc~^ gave 
the best consistency results for the stellar response to the spiral 
pattern jMartos et alJ 2004). But we will explore the range 15- 
30 km s~^kpc~^, assuming rigid rotation, as deduced from several 
literature determinations ( GerhardlboiOl) . This range gives a rota- 
tion period between 400 and 200 Myr. 

Density contrast. Determinations for the density contrast of 
the MW spiral structure are few and entail a large uncertainty. 
Moreover, the density contrast definition is sometimes ambiguous. 
For external galaxies th e spiral amplitude is often quantified as 
A2 (lGrosb0l etal.ll2OO4l) . which is the amplitude of the m — 2 
component of the Fourier decomposition of the surface brightness 
scaled to azimuthally averaged surface brightness. A different mea- 
sure is the arm-interarm c ontrast K{R) = (1 + A2)/(l — A2) 
teinnev & Tremain32008h . If A2 is measured fro m infrared b ands, 
a mass to light ratio of ~ 1 can be assumed jKentl [l993) and 
these quantities are directly a measure of the density contrast A2 ^ 
5a /ao (or K — (ctq + 5(t)/{(Tq — da)) where ctq is the axisym- 
metric surface density, and ^cr is the enhancement of density on the 
spiral arms. 

Table [2] summaris es several determinations of the density con- 
trast in the literature. Drimmel & Spergel (2001) give an arm- 
interarm ratio in the K band surface brightness in our Galaxy of 
K = 1.32 (or, equivalently, A2 = 0.14). However, these authors 
report that this may be undervalued due to underestimation of the 
arm scale height as compared to that of the disc. Indeed, this con- 
trast is significantly smaller than in external galaxies, where values 
up to A 2 — 0.6 are found ( Rix & ZaritskvJ995) . Recent data from 
the GLIMPSE survey give an excess of 20-30% of stellar counts at 
the maximum, with respect to the axisymmetric exponential fitting 
(Benjamin et al. 2005). For the exponential fitting, these authors 
have subtracted the enhancements in the spiral arms. Therefore, we 
ascertain that their value is a ratio with respect to the minimum 
stellar counts. Although the conversion to density contrast is not 
straightforward, we will assume that this gives directly K ^ 1.3 
{A2 ^ 0.13), which is not far from the Drimmel & Spergell (Eool 
value. We will take = 0.14 as a lower limit for the MW from 
iDrimmel & SpergellfcOOlh . 

Another determination for the MW spiral density contrast can 
be deduced from the relation between pitch angle and spiral ampli- 



4 Antoja et al. 



tude for external galaxies, explored in lGrosb0l et al.l (l2004h . For a 
pitch angle of 15.5°, the contrast A2 ranged approximately from 
0.1 and 0.23. Galaxies with a contrast up to 0.5 were also found, 
but the authors explain that this may be overestimated due to the 
strong star formation in these galaxies. Therefore, we will assume 
an upper limit of the density contrast of A2 — 0.23. 

Instead of being measured globally or as a function of radius 
as in external galaxies, for our own Galaxy density contrast esti- 
mations come from data of certain positions of the disc (mainly 
spiral tangencies, i.e. at a given radius). In the absence of detailed 
information about the density contrast along the MW arms, we as- 
sume that the chosen density contrast range is for an intermediate 
radius of ^ 6 k pc. Besides, we adopt an exponential fall as in 
IContopoulos & Grosb0l (1986.) with scale length of i?E = 2.5 kpc. 



3 MODELLING THE SPIRAL ARM GRAVITATIONAL 
POTENTIAL: TWO DIFFERENT APPROACHES 

Now we deal with the models for the gravitational potential of 
the arms. One approximation to model spiral arms would be the 
one given by Gerola & Seiden (197 5), who suggested that self- 
propagating star formation in a differentially rotating disc is ca- 
pable of producing large scale spiral features. This description pro- 
duces stochastic spiral arms and, although this may be able to ex- 
plain flocculent spiral arms, it would neither be the case for about 
the hal f of galaxies that pre sent grand design spirals in the infrared 
bands jKendall et alj|201l]) nor the case of the MW where clear K- 
band arms are observed (Section |2). Therefore, we do not consider 
a gravitational model for the flocculent arm type. 

The first model that we consider in this study is the TWA spi- 
ral arms. This density wave theory (Lin, Yuan & Shu 1969) has ex- 
tensively been used in stellar and gas dynamics simulations. From 
its birth, this spiral structure theory has been capable to provide 
certain qualitative explanation, going from some local kinematics 
of our own Galaxy to systematic changes in spiral arm properties 
along the Hubble sequence. However, it has never been straightfor- 
ward to obtain a good fitting between density wave theory, and ob- 
servations, especially, of course, when galaxies reach a non-linear 
regime, which might be the case in the majority of galaxies. One 
simple example among several: the density wave theory says that 
stronger shocks of material (gas, stars) with spiral arms should be 
observed in thinner spir al arms, but it seems th at wider arms pro- 
duce stronger shocks (Kennicutt & Hodgel ll982'). 

The spiral arms obtained through the TWA are an elegant so- 
lution to stability analysis of gal actic discs, by p erturbing the ba- 
sic equations of stellar dynamics dLin et al.ll 19691) . Treated as weak 
perturbations (low mass, low pitch angles) to the background po- 
tential, a simple cosine expression for the s piral arm potent ial is a 
self-consistent solution for the linear regime dLin et al.l 19691) . Nev- 
ertheless, it has been a costume to extrapolate this mathematical 
approximation from the linear domain to all kind of spiral galax- 
ies, under the general assumption that spiral arms are weak pertur- 
bations in galaxies. A key question then is whether the TWA can 
be applied to the MW spiral arms. First, the TWA solution is lim- 
it ed to tightly wound spiral a rms with m/ tani >> 1. As stated 
in lBinnev & Tremaine (2008), this is satisfied in most galaxies but 
not with a comfortable margin. According to the observational con- 
straints for our Galaxy (Section|2), pitch angles of i = 15.5° (locus 
1) and 12.8° (locus 2) give 7.2 ^ m/tani ^ 8.8, and, therefore, 
the assumption is at least doubtfully satisfied. Second, the TWA 
is a perturbative solution for small density contrasts. Although the 



determinations for the MW spiral density contrast entail large un- 
certainty, the maximum spiral arm density contrast of 23% of the 
axisymmetric disc seems to exceed the requiment. Due to the un- 
certainity on the MW spiral arm characteristics, it is not clear yet 
whether they can be compatible with the self-consistent TWA solu- 
tion. 

The second studied model for the spiral arm gravitational po- 
tential is the PERLAS model (sPiral arms potEntial foRmed by 
obLAte Spheroids, Pichardo et al. 2003). Unlike the cosine poten- 
tial of the TWA, which represents a local approximation to the 
spiral arms, the nature of PERLAS is a very different one. It is 
a model with material arms in the sense that it corresponds to a 
given spiral arm mass distribution, from which the potential and 
forces are derived. In particular, the model is constructed as a su- 
perposition of small pieces of mass distribution. Because of this, 
it is more flexible than the TWA. It can be easily adjusted to the 
MW spiral arms or any spiral density profile that might be far 
from following a cosine function , as the profiles in some exter- 
nal galaxies (K endall et alj|20li]) . For this model, we can adjust 
the total spiral mass, the arm width and arm height, which is not 
possible for the TWA model. This potential results more realistic 
in the sense that it considers the force exerted by the entire spiral 
arm structure, sculpting much more complicated shapes in poten- 
tial and force than a simple cosine function. These intrinsic differ- 
ences may induce significant deviations on orbital dynamics from 
the classic cosine. In addition, it is a full 3D model which, instead 
of taking an ad hoc dependence on the z coordinate (e.g an ad- 
ditional potential term sech^[z/^s]), it considers directly a three- 
dimensional mass distributio n. It also satisfies a p eriodic orbit di- 
agnostic for self-consistency dPichardo et al.ll2003l) , which consists 
in analysing the stellar orbital reinforcement of the potential as in 
Patsis, Contopoulos & Grosb0l ( 1991). As for the TWA, the orbital 
self-consistence of the model is only assured for a certain range 
of parameters, especially mass and pattern speed. Next we present 
in detail these two potentials and how they are tuned to the recent 
observations of the MW spiral structure. 



4 THE MODELS 

The axisymmetri c part of our MW models is taken from 
lAUen & Santininl (Il99lh (hereafter A&S). It is composed by a 
bulge, a flattened disc and a massive spherica l halo. The first two 
are m odelled as Miyamoto-Nagai potentials jMivamoto & Nagall 
Il975h and the halo is built as a spherical potential. The main 
adopted observational constraints of the model are summarised in 
table 1 of A&S. A value of Rq = 8.5 kpc for the Sun's galacto- 
centric distance and a circular speed of ^(^^0) = 220 km are 
adopted. The total axisymmetric mass is Mr = 9 x 10^^ M©. The 
local circular frequency is Q.{Rq) = 25.8 kms~^kpc~^. The two 
spiral models, TWA and PERLAS, are described in next sections. 

4.1 The TWA spiral arms 

In the TWA the spiral arm potential is obtained by solving 
the Poisson's equation for a plane wa ve in a razor-thin disc 
dBinnev & Tremainel2008l : lLin et al.ll969h . after modelling the spi- 
ral structure as a periodic perturbation term to the axisymmetric 
disc. Following IContopoulos & Grosb0ll (Il986 ), the TWA gives a 
potential in the plane of the form: 

$,^(i?, 0) = -A,pi?e-^/^^ cos(m[0 - 0o] + g{R)), (2) 
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Figure 2. Density contrast A2 = Sa/cro as a function of radius of the TWA 
model for the two ampHtudes that fit the high and low observational limits 
(vertical error bar) for the MW. 
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Figure 3. Density contrast K = (ao + 6a) /(ao) — 5(7 as a function of 
radius of the PERL AS model for the two amplitudes that fit the high and 
low observational limits (vertical error bar) for the MW. 



Table 3. Assumed parameters of the TWA model. 



Table 4. Assumed parameters of the PERLAS model. 



Parameter 






Range 


Amplitude (locus 1) 


Asp {[kms ^ 


2kpC-l) 


850-1300 


Amplitude (locus 2) 


Asp ([kms~^ 


2kpC-l) 


650-1100 



where Asp is a measure of the amplitude of the spiral pattern, 
and the rest have already been defined. T he amplitude of the spi- 
perturbing potential, following again IContopoulos & Grosb0ll 
(1986), it is related to the amplitude of the perturbed surface density 
(through Poisson's equation) as: 



6a{R) = 



-R/Rj: 



7vG\ tani| 



(3) 



To fit this model to the MW spiral arms we use the locus and 
parameters as discussed in Section |2] In particular for this model, 
the density contrast A2 = 6a(R)/ao(R) is obtained directly from 
Equation ^ and the surface density ao of the stellar part of the 
axisymmetric A&S model (disc and bulge). Due to the different 
dependence with R of 6a and the disc surface density in A&S, the 
density contrast A2 decreases with radius. As explained in Section 
[21 we chose the range of Asp that fits the observational limits for the 
density contrast A2 at a intermediate radius ofR^G kpc. In Table 
[3] we show the determined range. For locus 1, this is Asp =850- 
1300 [ km s~^]^ kpc~^ . The density contrast as a function of radius 
for these two limits is shown in Fig. |2] For locus 2 (with a smaller 
pitch angle) the range is Asp =650-1100 [kms-^]^kpc-\ 

Notice also that Equation ([T]) gives a locus that joins the two 
starting points of the spiral locus in the disc central part (not plotted 
in Fig.[T]). However, we have checked that our simulations (which 
are studied in more outer regions) do not depend on the exact shape 
in these inner parts. 



4.2 The PERLAS spiral arms 

The spiral arms of the PERLAS model consist of a mass distri- 
bution which is built as a superposition of inhomogeneous oblate 
spheroids along a given locus. A linear density fall, inside each 
spheroid, is considered with zero density at its boundary. The cen- 
tral density of the spheroids follow an exponential fall in R along 
the arms with scale length of T able [T] The po tential and force fields 
for these spheroids are given in ISchmidtl([T956i) . The overlapping of 
spheroids allows a smooth distribution along the locus, resulting in 
a continuous function for the gravitational force. We checked that 



Parameter 



Values or range 



Beginning of the arms 
End of the arms 
Arm half- width 
Arm height 
Mass (locus 1) 
Mass (locus 2) 



Ri (kpc) 
Rf (kpc) 
ao ( kpc) 
CO ( kpc) 
Msp/Md 
Msp/Md 



3.3 

12 
1 

0.5 
0.03-0.05 
0.035-0.06 



no significant change was observed if this separation was decreased 
thus increasing the smoothness of the spiral mass distribution . For 
more details about the construction of the model see Pich ardo et al.l 
(l2003h . 

To fit this model to the MW spiral arms we use the locus and 
parameters as discussed in Section|2] Additional parameters for this 
model are shown in Table |4] In this case, the spiral amplitude is 
quantified through the value Msp/Mo which is the ratio of the spi- 
ral mass to the mass of the disc in the A&S model. In the PERLAS 
model the spiral arms are added as an mass enhancement to the ax- 
isymmetric background on the imposed locus Q Because of this, for 
the PERLAS model the mass of the spiral arms {Msp) is globally 
subtracted from the original disc of the A&S model (Md), which 
guarantees that the total mass of the model do not change and nei- 
ther does the mean circular velocity. For this reason, the parameter 
K is more suitable to be related with the Msp/Mn, as it is the ra- 
tio of the surface density on the spiral arm (axisymmetric disc and 
spheroids) to the minimum density (axisymmetric disc). As for the 
TWA model, the density contrast K decreases with radius. For lo- 
cus 1, the range of spiral mass ratio that fits the observational range 
of density contrast K Sii R ^ 6 kpc is Msp/Mn =0.03-0.05 (Fig. 
[3]). For locus 2, the respective range ifl Msp/Mo =0.035-0.06 

Additionally, the beginning of the spheroid superposition, i.e. 
the effective arm beginning, is at 3.3 kpc and 3.6 kpc, respectively 
for locus 1 and locus 2. The superposition ends at 12 kpc for both 

^ Note the difference with the TWA for which the arms consist of a density 
perturbation that adds density in the arm region but subtracts density in the 
interarm region. 

^ Contrary to the the TWA model where a smaller pitch angle demanded 
a lower amplitude to reproduce a given density contrast, for the PERLAS 
model a higher mass ratio is needed. This is because we fix the spiral arm 
end at = 12 kpc and for locus 2 with a more tightly wound arm, the arm 
longitude is larger than for locus 1 . This implies that a higher total spiral 
mass is necessary to reproduce the desired contrast. 
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Figure 4. Decomposition of the radial force of the PERLAS model in the 
terms Qmr for m = 2, 4, 6, 8 of Equation f^. Long dashed line is the ratio 
of the radial axisymmetric part of the force of PERLAS model to the total 
axisymmetric radial force <I>q gp/^Q. 



cases. The arm half- width and height above the plane are taken 
to be 1 kpc and 0.5 kpc, respectively, according to the analysis of 
a sample of external galaxies (Kennicutt & Hodge. 1 982.) and the 
discussion in lMartos & Coxl (Il998h . 



5 FORCE FIELDS OF THE SPIRAL MODELS 

The main difference between the PERLAS spiral arms and the 
TWA falls on the construction itself. In the TWA, the spiral arms 
are a small perturbative term of the potential. This produces a can- 
cellation of the contribution from the distant parts of the pattern to 
the local force. The PERLAS spiral arms correspond to an inde- 
pendent mass distribution. The contribution from the entire spiral 
pattern causes the spiral potential and force to adopt shapes that are 
not correctly fit by the simple TWA perturbing term that has been 
traditionally employed. In this section we show these fundamental 
differences in detail by studying the force field exerted by these two 
models using the Fourier decomposition. If not stated the contrary, 
the plots and values in this section are referred to the lower limit of 
the density contrast defined in Section |2] 

First we study the parameter qmr, which is the m term of 
Fourier decomposition of the non-axisymmetric radial force scaled 
to the axisymmetric radial force. This is: 



qmr{R) — 



(4) 



where primes denote derivatives with respect to R, and ^rn{R) 
is the m-component of the Fourier decomposition of the global 
potential, that is of the potential of the spiral arms plus the A&S 
potential. In particular for m = 0, $o (R) corresponds to the ax- 
isymmetric part of the global potential. The TWA model consists 
of a pure and simple m = 2 component. Therefore, the spiral arms 
only contribute in the numerator of Equation (|4), and the axisym- 
metric model A&S is the denominator. By contrast, the PERLAS 
arms consist of a more complex potential structure with more than 
a m = 2 term. In Fig. |4] we show the Fourier decomposition of 
the PERLAS potential. In particular, the long dashed line shows 
the ratio of the axisymmetric component of the radial force of the 
PERLAS model ^o^sp to the axisymmetric global radial force $o. 
We see that the axisymmetric radial force of the PERLAS model 
can be as large as 4% of the total axisymmetric force (and up to 
7% for high limit of the density contrast). This important m = 




Figure 5. Radial spiral force qr as a function of radius for the TWA and the 
PERLAS model. 



component contributes, toguether with the A&S model, to the ax- 
isymmetric part of the global potential (denominator in Equation 
(lU) together with the A&S model. In the same plot we see that 
for this model the m = 2 component is dominant but it also has 
non- vanished m > 2 terms qmr- 

In order to quantify the whole non-axisymmetric force of our 
spiral arm models we define the parameter: 



qr{R) - ^ qmr{R) 



(5) 



which is the ratio of all the non-axisymmetric radial force terms 
to the axisymmetric part of the global potentially For the TWA, 
q^ — q2r. But for the PERLAS model it is the contribution of all 
terms qmr, with m ^ 2, of Fig.|4] Fig. |5] shows qr as a function of 
radius for the two models. The range of the radial force amplitude 
of the two models is approximately the same as they are fitted to 
reproduce the same density contrast. However, the detailed shape 
of qr{R) is different due to the important contribution of all terms 
m ^ 2 in the PERLAS model. From this figure we also observe 
that the maximum radial force of PERLAS spiral arms with re- 
spect to the total axisymmetric radial force is 0.09 (absolute value). 
Around Rq the parameter qr is 0.003. For the TWA, the maximum 
is 0.14 and the value at solar radius is 0.03. For the higher limit of 
the density contrast range, the shape of qr{R) does not change. But 
in that case, the PERLAS model produces a maximum radial force 
with respect to the axisymmetric background of 0.16 and a value 
near the solar radius of 0.005. For the TWA the maximum value is 

0.21 and the value at the solar radius is 0.04. 

The parameter Qt{R) was used in 'Combes & Sandersl(ll98lh 
to quantif y the tangen tial force of bars. Here we use it for spiral 
arms as in lBlock et al It is the ratio of the maximum spi- 

ral arm tangential force at a given radius to the mean total radial 
axisymmetric force at that radius: 

The parameter Qt includes the tangential forces due to all terms 
m ^ 2 for the PERLAS model. This parameter is shown in Fig. 
[6] The large and small bumps at ^ 3.5 kpc and 12 kpc for 
the black curve are due to the more abrupt beginning and end of 
the arms in the PERLAS model but have no consequences for our 



^ Note that often, e.g. in lAthanassoula et al] (ll983V the strength of the non- 
axisymmetric components is quantified through the parameter qr which, 
contrary to the present study, only includes the m = 2 component. 
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Figure 6. Maximum tangential spiral force Qt as a function of radius for 
the TWA and the PERLAS model. 
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Figure 7. Radial force (top) and tangential force (bottom) as a function of 
radius for two different azimuths for the TWA and PERLAS model. The 
forces are scaled to the A&S radial axisymmetric force. 



Study that is focused on near solar radius. The amplitude range of 
Qt is rather similar in the radii of interest in this study (7-10 kpc). 
The maximum of Qt{R) gives a single an d quantitative mea sure 
of the torque or strength of the spiral arms (block et al.ll2004h and 
is called Qs. For the PERLAS model we find that = 0.1 and 
for the TWA it is 0.039. For the maximum density co ntrast, Qs is 
0.17 fo r the PERLAS model and 0.060 for the TWA. In lBlock et aP 
(I2OO4I) the spiral strengths Qs in a sample of 15 external galaxies 
range from to 0.46. According to this, the models for the MW spi- 
ral arms are below the median of external galaxies. For both models 
Qt{Rq ) is 0-01 for the lower density contrast limit and ^ 0.02 
for the high case. 

The parameters qr and Qt measure average and maximum 
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values at a given radius, respectively. We now study the detailed 
force profile. Figs. |7] and [8] show the radial and tangential forces of 
the two models as a function of azimuth and radius R scaled to 
the axisymmetric radial force of the A&S model. We see signifi- 
cantly different force fields. In general, the PERLAS spiral arms 
(solid line) presents more abrupt features than the TWA. It also 
has different locations (R, (j)) for the minima and maxima for both 
the radial and tangential forces. The positions at which the force 
changes its sign is different as well. Besides, the radial force of the 
PERLAS model is in general below the force for the TWA, that 
is non- symmetric with respect to and shifted to negative values. 
This is due to the inner enclosed spiral mass at each radius, that 
is the m = component of the PERLAS spiral force mentioned 
before. For instance, in the top right panel in Fig. |8] (solar radius 
i? = 8.5 kpc) the radial force is negative for all 0, whereas the 
TWA consists of a symmetric oscillation around force. 

To conclude, it is quite different to model the gravitational po- 
tential by using the TWA (i.e. a cosine function, that fits a given 
locus and certain density contrast) or by building a mass distribu- 
tion on the same locus and contrast. We have seen that the two 
approaches to model the spiral arms gravitational potential give 
significantly different force profile. In next sections we study how 
these differences propagate to the kinematic stellar response. 



Figure 8. Same as Fig. |7] but as a function of azimuth for two different 
radius R. 



6 THE SIMULATIONS 

In order to study the effect of the spiral arms on the kinematic dis- 
tribution near the Sun, we perform numerical integrations of test 
particle orbits, as most of the studies up to now. These are sim- 
ple models and their self-consistency is not assured (Section[3]). By 
contrast, models such as N-body simulations of galaxy formation 
could model self-consistently the spiral arm kinematic effects, as 
well as include naturally the bar or spiral arm evolution and other 
processes that might also sculpt the velocity distribution (e.g. past 
accretion events, star formation bursts). However, N-body simula- 
tions with larger than current number of disc particles, better spa- 
tial and temporal resolution, and models similar to the MW are 
required. At this moment, test particle simulations allow us to use 
models that fit the MW spiral arms in a controlled manner. They 
are simpler and faster models that offer easy exploration and un- 
derstanding of the influence of the spiral arm properties. 

For our simulations we adopt the potential models of Section 
IHand the initial conditions described in Section [6Jl In Section [Ol 
we give more details about the method. Our study focuses on the 
kinematic effects produced on or near the plane. For this we assume 
that the vertical motion is decoupled with the in-plane movement. 
This is reasonably true for ne arly circular orbits that do not take 
larger height above the plane (iBinnev & Tremainell2008l) . As our 
initial conditions consist of rather cold discs (see Section [SjI , we 
adopt these assumptions and simplify our analysis by considering 
2D simulations {z — 0). 

The ([/, V) velocity reference system is used which is centred 
on a given position on the disc plane and moves following the Re- 
gional Standard of Rest (RSR). This is defined as the point located 
at a radius R that describes a circular orbit around the Galactic 
Centre with a constant circular speed Vc{R). U is the radial veloc- 
ity component, which is positive towards the Galactic Centre, and 
V is the tangential component, positive in the direction of Galactic 
rotation. 
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6.1 The initial conditions 

We explore two different types of initial conditions: ICl & IC2. 
Both are axisymmetric discs, following an exponential profile with 
scale length of i?s = 2.5 kpc. This is similar to the thin disc scale 
length found by Freudenreich (1998) from COBE/DIRBE data. 
ICl- The velocity distribution relative to the RSR is adopted as a 
Gaussian with low dispersions au — cry — 5kms~^ and constant 
for all radii. These values are similar to the induced gas veloc ity dis- 
persions in the plane due to Galactic spiral shocks (Kim 2009) and 
to the dispersion of the youngest Hipparcos stars ( Aumer & Binnev 
12009). With this IC we aim to simulate a cold young disc. 

IC2- For these IC, the p hase space distribu tion function (DF) 
is constructed, as discussed in lHemQuistl(ll993h . The velocity field 
is approximated using the moments of the collisionless Boltz- 
ma nn equation, sim plified by the epicyclic approximation. Follow- 
ing |DeMe3 (ll999D , we will adopt proportional to the surface 
density, and therefore, the radial velocity dispersion profile is: 

at7(i?) = aoe-^/('^^\ (7) 

The local normalisation is chosen au (Rq) ^ 20 km . The tan- 
gential velocity dispersion profil e is derived from the epicyc lic ap- 
proximation (equation 4.317 of iBinnev & Tremainel l2008h . Also 
we take into account the asymmetric drift (equation 4.228 of 
iBinnev & Tre maine 2008). Assuming that the disc is stationary, ax- 
isymmetric and symmetric about its equator, and that the orienta- 
tion of the velocity ellipsoid is aligned with the coordinate axes, it 
is: 

"• = l|[l^-' + «(i + 5i)l- '« 

This DF confers to the IC2 disc the properties of an in- 
termediate population (^ 1 Gyr) o f the MW thin disc 
(iHolmberg. Nordstrom & Andersen[l2009h . 



adjacent U-V planes. The symmetries of the Galactic potential 
($(i?, (/))=$ (i?, 180° + 0)) allow us to double the number of par- 
ticles in each studied region. In all cases, the number of particles in 
the final velocity distributions is larger than ^ 10000, being statis- 
tically robust. Fixed the spiral arm pattern speed, these 15 regions 
have different ratio Qsp/^, that is, the pattern speed scaled to the 
local azimuthal frequency of that region. These are representative 
of different pattern speeds, allowing us to explore small variations 
in this parameter. Also they have different relative spiral phase (l)sp 
which aims to explore the range of uncertainty in this parameter, as 
explained in Section |2l 

We apply the Wavelet Denoising method (WD) to the velocity 
distributions. This method gives a smooth DF from a discrete point 
distribution via a smoothing/filtering treatment at dif ferent scales, 
that e liminates Poisson fluctuations (for details, see lAntoia et"aD 
l2008h . The results will be shown with logarithmic colour scale and 
contours representing, from inside out, 0.01, 0.05, 0.10, 0.20, 0.30, 
0.40, 0.50, 0.6, 0.70, 0.8, 0.9 and 0.999 of the maximum density. 



7 STELLAR RESPONSE TO THE TWA AND PERLAS 
MODELS 

We have run simulations for pattern speeds Qsp of 15, 18, 20, 22, 25 
and 30 kms~^kpc~^ (inside the observational range of evidence, 
Section|2]) for both models. We have used the two limits of density 
contrast, the two loci and the different initial conditions detailed in 
previous sections. Here we compare the results for the two models. 
First, in Section ITJI we analyse the density response to the spiral 
models in our simulations. With this we aim to study the consis- 
tency between density response and gravitational potential of the 
models, and to relate the density response to the kinematic analysis 
of subsequent sections. 



6.2 Integration procedure and analysis 

In our simulations the integration of each particle is initialised at 
a time t = — r (as a convention) and ends at t= 0, being r the 
particle exposure time to the potential. The time r is chosen at ran- 
dom between and 2 Gyr for each particle. This maximum inte- 
gration time corresponds to ^ 7 revolutions of the spiral arms for a 
typical pattern speed of 20 km s~^kpc~^ . With this procedure, the 
final velocity distributions result of a superposition of particles in- 
tegrated different times, resembling the observed distribution, with 
a superposition of stars of different ages. In some cases, we sep- 
arate the particles into different bins of integration times, to study 
the induced kinematic effects as a function of time. 

The integration of the motion equations is done with the 
Bulirsch-Stoer algorithm of IPress et al.l (1 19921) . conserving Ja- 
cobi's integral within a relative variation of | {Eji — Ejf) / Eji \ ^ 
10~^. The reference frame used for the calculations is the rotation 
frame of the spiral arms. The number of particles in each simulation 
is about 10^. 

After the integration, we study the induced kinematic distri- 
bution near the Sun. For the analysis, we focus on the U-V plane 
of the 15 circular regions indicated in Fig. [T] which have 300 pc 
of radius (similar to the radius of the available observed local ve- 
locity distribution, see fig. 1 in Antoja et al. 2008). The centres 
of the regions are located at five different radius (7.3, 7.9, 8.5, 
9.1 and 9.7 kpc) separated 600 pc from the nearest ones, and at 
azimuths —15°, 0° and 15°, which is found to be the approxi- 
mate azimuth interval that shows significant differences between 



7.1 Density response 

Fig. |9] shows the density distribution of a simulation with the TWA 
(top) and the PERLAS model (bottom) taking the same conditions 
(pattern speed, locus, etc.) for both models. Darker regions show 
the denser regions with the same scale in all panels of the figure. 
First and second column correspond to short and long integration 
time, respectively. Blue stars show the positions of maximum den- 
sity in rings of 500 pc of radius. The white curves indicate the locus 
of spiral arm model. Inner and outer red dashed circles show the 
approximate positions of the 4:1 inner resonances and corotation, 
respectively. 

First row shows that the density response to the TWA with 
the minimum density contrast and a pitch angle of i = 15.5° does 
not follow the spiral locus imposed in the potential approximately 
beyond the 4:1 inner resonance. The spiral arms in this simulation 
show a more tightly wound pattern at outer radius for short inte- 
gration time and high dispersion for long integration times. Right 
panel also shows complex structures like rings and low density re- 
gions. The different pattern speeds and amplitudes in the ranges 
esta blished for the MW give equival ent results. 

IContopoulos & Grosb0ll (Il986h already studied the density re- 
sponse for several spiral models based on the TWA, mainly focus- 
ing on periodic orbits. They concluded that strong spirals can not 
exist beyond the 4:1 inner resonance as in this disc region the spi- 
ral response is out of phase due to non-linear effects. By contrast, 
weak or tightly wou nd spirals can e xist up to corotation, as the lin- 
ear theory predicts (iLin et al.l 1 19691) . Our results are in agreement 
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Figure 9. Comparison between the imposed spiral locus and the density response to the TWA model (top) and PERLAS model (bottom) for Qsp = 
18 km s~^kpc~^, locus 1, ICl and minimum density contrast. The density is scaled to the axisymmetric exponential density at the beginning of the simula- 
tion. First column shows the simulation for small integration times (0-400 Myr) and second column corresponds to large integration times (1600-2000 Myr). 
Red dashed circles show the approximate positions of the 4:1 inner resonances (inner circle) and corotation (outer circle), estimated from the axisymmetric 
model of A&S. 



with the results of 'Contopoulos & Grosb0ll ('1986') for strong spi- 
An exhaustive analysis of the response density is out of the 
scope of this paper. But the difficulties in using the TWA tuned to 
the MW as a self-consistent model beyond the 4: 1 inner resonance, 
poses important constraints in the use of the TWA for the MW. For 
certain pattern speeds, this would imply a too short ending radius 
of the spiral arms (e.g. the spiral arms could end at ^ 4.3 kpc for a 
pattern speed of 30 kms~^kpc~^)|3To conclude, the TWA might 
not be a proper description for the MW spiral arms for some pa- 
rameter combinations (pattern speed, pitch angle, amplitude). 

The density response to the PERLAS model (second row of 
Fig. 17. Il l shows also disrupting effects near the inner 4:1 resonance 

^ Preliminary tests seem to indicate that the non-linear effects in our simu- 
lations are due mainly to MW spiral arms that are not enough tightly wound. 
A simulation with the same pitch angle, but smaller amplitude, still presents 
the inconsistency between imposed pattern and response. On the contrary, a 
simulation with the same amplitude but smaller pitch angle follows almost 
perfectly the locus up to corotation for all integration times. 
^ Notice that the existence of an additional slower spiral mode at outer ra- 
dius could solve this problem. We are not considering this type of composite 
models in the present paper, but see discussion in Section [TL2] 



but, at larger radius the density response does follow the imposed 
spiral arms for small and large integration times. Indeed, with larger 
pattern speeds, the response follow the imposed locus even be- 
yond corotation0 Nevertheless, the orbital self-consistence of the 
PERLAS model is only assured for a certain range of parameters, 
specially in mass and pattern speed (see Pichar do et al.l (1200 3l) ). 

For the PERLAS model we also see complex density struc- 
tures such as two extra arms for short integration times or low den- 
sity regions and rings in the long integration case. All these patterns 
that are observed for both models could be related to the different 
resonances of the spiral arms and might well represent rings and 
spiral arm spurs that are seen in our Galaxy and in external galaxies 
jElmegreen & Elmegreen|[l99Ql) . For instance, two emptier regions 
are located on the 4:1 inner resonance at i? ^ 8.0 kpc in the right 
panels. Notice, however, that these appear at different azimuths de- 
pending on the model (near = 0° for the PERLAS model but at 



^ It i s worth mentioning that other models different from the TWA (iToomrel 
"1981"; 'Rautiainen & Sald'l999|; IVogUs. Tsoutsis & E fthvmiopo ulosl 1200^ 
Romero-Gomez et al. 2007) have shown that spiral arms can extend well 
up to or beyond corotation as it happens with the PERLAS model. 
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Figure 10. Zoom of the density response at solar position for the TWA 
model (top) and PERLAS model (bottom) for all particles (all integration 
times) in the simulation of Fig.|9] 



an azimuth of ~ 50° for the TWA). In fact, note that the induced 
density in our region of study (near solar positions) is significantly 
different between models. Fig.[TO]is the local density response for 
all particles of the simulation. The TWA model (top) presents softer 
density gradients, whereas the PERLAS model (bottom) shows the 
mentioned emptier region at solar azimuth. The local density re- 
sponse of the two models is different in all our examined cases, 
emphasizing the distinct force fields studied in Section |5] and fore- 
casting the differences in the induced velocity distribution near so- 
lar positions (Section 1X21 ). 

Simulations with IC2 show qualitatively equivalent results to 
the ones with ICl for both models, except for the fact that the den- 
sity features appear less sharp but more diluted, as expected for 
a hotter and, therefore, less responsive disc. For both models we 
observe a transient nature of the density patterns. In all cases the 
density response becomes stationary approximately between 400 
Myr and 1200 Myr of integration time depending e.g. on the pat- 
tern speed (except for outer radius R > 10 kpc, out of the region of 
study). Similar time dependence will be observed for the kinematic 
distributions of these models (Section [TO]). 



7.2 Different kinematic imprints 

In this section we compare the stellar kinematic response to the 
PERLAS model and the TWA. As an example, we now adopt 
locus 1, maximum density contrast, ICl and a pattern speed of 
18 kms~^kpc~^ for both cases. The U-V plane for the regions 
of Fig. [Dare shown in Figs. [TT] and [12] for the TWA and PERLAS 
models, respectively. The 15 panels of each figure are positioned 
similarly to their spatial location of Fig. [T] Each panel corresponds 
to a region with different relative spiral phase (/)sp and ratio Qsp/^, 
which are indicated in each panel. 

First we see that both models induce several and rich kine- 
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Figure 11. U-V velocity distributions for the 15 regions of Fig.[T]for the 
simulations with the TWA model with locus 1 , maximum density contrast, 
ICl and a pattern speed of 18 km s~-'^kpc~-'^ . The 15 panels are positioned 
in the figure similarly to their spatial location of Fig. [T] The scale in the 
right indicates the spatial position of the spiral resonances (see text). 



matic substructures near solar position. Nevertheless, we see that 
the details and induced groups at each of the 15 regions are dif- 
ferent for each model, in agreement with the differences between 
induced density (Section ITJI . For instance, an elongated group at 
low V appears for the left panel of the last row for the TWA model 
but it is present in all panels in this row (several spiral phases) for 
the PERLAS model. The central part of these panels is also fairly 
different. The rest of the panels also show groups at different posi- 
tions depending on the model. 

We have detected that the two models occasionally give a sim- 
ilar U-V plane when comparing different disc positions with com- 
parable density distribution. For instance, we see similar velocity 
distributions near the emptier region on the 4:1 inner resonance 
(Section [TTJ, which is is located at different azimuths depending 
on the model. We attribute this occasional shift to the differences 
in the force fields of both models (Section [5]). We have seen that 
maxima, minima and points of the force profiles (both tangential 
and radial) are locally shifted in radius and azimuth when compar- 
ing the two models. For example, in Fig. [8] we see that the max- 
imum tangential force of the PERLAS model is shifted to larger 
azimuths compared to the TWA. On the contrary, the minima are 
shifted to smaller azimuths. The link between phase space, force 
profile and resonance locations involves a large-scale study of the 
models which is out of the purpose of the present paper, but is cur- 
rently under study (Antoja et al., in preparation). 

As common aspects of the kinematic response of the two mod- 
els, we notice variety of shapes, sizes, inclinations and positions of 
the induced substructures, for this and other pattern speed ratios. 
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Figure 12. U-V velocity distributions for the 15 regions of Fig.[T]for the 
simulations with the PERLAS model with locus 1, maximum density con- 
trast, ICl and a pattern speed of 18 km s~-'^kpc~-'^ . 



For example, notice more or less rounded groups but also thin elon- 
gated structures. These and other simulations have shown that both 
models induce structures at V as high as 30 kms~^ and as low as 
— 60kms~^. The U component of the induced kinematic planes 
can range from -80 to 60 kms~^. The number of induced groups 
depends on the region but is in most cases 2 or 3. Interestingly, the 
kinematic groups are found to be in general more symmetric with 
respect to the [/ = axis for the TWA than in the PERLAS model. 

More important, by inspecting these kinematic plots, we see 
that the PERLAS model induces kinematic substructure where the 
TWA model does not (e.g. first and second rows). The TWA gives 
substructure for a smaller range of pattern speeds. We examine this 
in detail in Section [9l In the following sections, we study the influ- 
ence of the spiral arm properties on the kinematic distribution near 
the solar position. We focus basically on the PERLAS model but 
we contrast the results with the ones obtained by the TWA model 
when important differences arise. 



8 EFFECTS OF THE PATTERN SPEED AND SPIRAL 
ARM ORIENTATION 

Figs.[T3land[T4lshow the contour plots of the U-V plane for the 15 
regions shown in Fig.[T]of the simulations corresponding the pattern 
speeds of 15 kms~^kpc~^ and 20 km s~^kpc~^, respectively, for 
the PERLAS model. Together with the already presented Fig. [12] 
these panels allows us to study the effects of several pattern speed 
ratios Q^sp/^ and relative spiral phases (j)sp. In all cases we use 
locus 1, maximum density contrast and ICl. 

The panels show rich kinematic substructure near the solar 
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Figure 13. Same as 
15 km s^-'^kpc""'^. 
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Figure 14. Same as 
20 kms-^kpc-^ 
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Figure 15. Same as Fig.[T2lbut for locus 2. 



Figure 16. Same as Fig.[T2]but for the initial conditions IC2. 



position, strongly depending on the pattern speed ratio and relative 
phase of the spiral arm. The kinematic substructures change more 
slowly with azimuth (along rows) than with radius (along columns). 
The changes are significant if one moves only ^ 0.6 kpc in radius. 
But 2 kpc are needed in azimuth to detect important differences 
in the U-V plane. The kinematic groups slightly change position 
on the U-V plane with azimuth or radius. For instance, for a pattern 
speed of ^sp/^ 0.65 (fourth row in Fig. [12) the velocity distri- 
butions show basically two groups, elongated in the U direction, 
which appear in different positions in the U component depending 
on the azimuth. 

There are equivalent U-V plane plots in different simulations. 
For instance, top left panel of Fig. [13] and third panel of the fourth 
row of Fig. [T2I correspond to similar relative pattern speed ^sp/^ 
and spiral phase (/)sp. However, moving in radius is not strictly 
equivalent to changing the pattern speed. First the spiral strength 
or density contrast depends on radius (see Section |4). Second, al- 
though the initial velocity dispersion of particles in ICl is flat with 
radius, at the end of the simulation this increases for inner radius, 
as we have seen with a simulation only including the axisymmetric 
part of the model. One example is the group at ~ — 40 km 
in the right panel in the last row Fig. [12] which is not present in 
the equivalent panel of Fig. [13] (left panel of the second row). The 
smaller strength of the spiral at outer radius and/or the smaller ve- 
locity dispersion do not allow to crowd this structure for the simu- 
lation in Fig. [13] 

Now we compare the results with locus 1 and 2 (Fig.[T). Fig. 
[T5] shows the same case as in Fig.[T2]but using locus 2. The small 
difference between the pitch angles of these loci (i — 15.5° and 
i — 12.8°, respectively), produces essentially the same kinematic 
structures when we compare regions with the same ratio Q^sp/^ 
and relative spiral phase (/)sp. Few differences are observed such as 



the group at low V for left panel of the last row in Fig.[T5]that does 
not appear in the equivalent relative position for locus 1 (approx- 
imately middle panel of last row of Fig. [12]). These subtle differ- 
ences are due to the tightness of the spiral arms for locus 2. 

Most of the 15 regions considered here correspond to inter- 
arm regions. When we analyse the velocity distributions in differ- 
ent regions of the Galactic disc, we see that the induced kinematic 
substructure increases near and on the spiral arms (Antoja et al., in 
prep.). For locus 2, we find more substructure than for locus 1 for 
the TWA because for this locus some of the regions are closer to 
the spiral arms. By contrast, for the PERLAS model and the same 
pattern speeds, rich substructure is found for both locus. 

As expected the hotter population defined by the IC2 initial 
conditions (Section 16.11 ) does not respond so strongly to the spi- 
ral perturbation. The main groups are still observed as the only 
changes are in the way each structure is populated but not in the 
orbital structure. The signal of some fine details slightly fades or 
appears less defined for the hotter case. We show an example in 
Fig. [16] Notice that the bi-modality in the first panel of the first 
row in Fig.[T2]that does not show up for IC2. Other structures man- 
ifested as separated groups are now stickled with particles filling 
the gaps in between. 



9 STRONG KINEMATIC SUBSTRUCTURE AND 
RESONANCES 

Changing the pattern speed Q^sp of the spiral arms corresponds to 
changing the position of the spiral resonances with respect to the 
regions under study. These resonances can be, as shown in other 
studies (lDehnen|[200Ql : lOuillen & Minchevll2005l) . a major influ- 
ence on the orbital structure. A detailed orbital analysis would be 
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needed to determine the exact influence of each of the resonance 
in sculpting the velocity plane. Here, we examine the effects of the 
resonance proximity to the 15 regions of study. In the right part of 
e.g. Figs. [12] [13] and [14] we show a scale in galactocentric radius 
(configuration space) where we indicate the approximate radius of 
the main resonances. We also show in this scale the limits of the 
regions corresponding to the 3 adjacent velocity distributions. 

We determine that clear and rich kinematic substructure is 
seen near solar azimuth for the PERLAS model for pattern speeds 
from 0.5 to 0.75 times the angular rotation rate. With a local circular 
frequency of 25.8 km s~^kpc~^, this is for pattern speeds from 13 
to 19.5 kms~^kpc~^ at solar radius. This range of pattern speeds 
corresponds to being near the inner 3:1, 4:1 and 5:1 resonances. 

We also see that for pattern speed ratio of 0.76-0.85 {Q^sp ap- 
proximately between 19.5 and 22 km s~^kpc~^) we detect minor 
groups, that is, few and close to each other, or only small deforma- 
tions of a unique clump centred in the U-V plane (e.g. second and 
third rows in Fig. [14). In these cases, the regions are located espe- 
cially near the 6:1 inner resonance or high order m:l resonances 
(with m > 6). Therefore, we attribute weaker effects to these res- 
onances. For pattern speed ratio approximately between 0.86 and 
1.2 {Vlsp between 22 and 31 kms~^kpc~^) we observe no sub- 
structure at the considered azimuths (e.g. upper row in Fig. [T4l) . 
This is near corotation and high order resonances (inner and outer 
m:l resonances with m > 6), concluding that these resonances 
have no effect on the velocity plane at solar azimuths and for the 
maximum spiral arm density contrast. Q 

If we compare the results for the maximum density contrast 
with the lower spiral mass limit, the velocity distributions are simi- 
lar but we see in most cases less populated or/and smaller kinematic 
groups for the low density case. As an example we examine in Fig. 
[TTl the same case as in Fig. [12] but using the lower observational 
limit for the density contrast. Notice how the structures are essen- 
tially the same but the groups in the extremes of the U-V plane (far 
from the (0,0) km point) disappear for the weaker spiral arms. 
For instance, the groups at low V of last row of Fig. [12] are not 
populated in Fig. [T7] We observe this trend also for other pattern 
speeds. The range of pattern speeds where we find significant kine- 
matic substructure is reduced with respect to the maximum density 
contrast case (0.5-0.75) and it now approximately 0.58-0.75. Thus 
being near the 3:1 inner resonances has lower influence for weaker 
arms. 

By examining more simulations, we see that the TWA gives 
less substructure than the PERLAS model, given the same density 
contrast. It also gives clear and rich kinematic groups near solar 
azimuth for a smaller range of pattern speeds. Specifically, for the 
maximum density contrast, this range is 0.60-0.73 (pattern speed 
approximately from 15.5 to 18.8 kms~^kpc~^). This is for re- 
gions close to the 4: 1 and 5:1 inner resonances. For smaller pattern 
speeds, only minor substructures or small deformations of the ve- 
locity distributions are observed, contrary to the PERLAS model. 
For larger pattern speeds, the induced kinematic groups disappear 
abruptly for the TWA, but progressively for the PERLAS model. 
For a fixed pattern speed, all these lead to a smaller range of radius 



^ Note that important substructure is also observed below the ratio 0.5 but 
is not considered as it is out of the established observational range. We 
also see that some substructure or deformation of a unique group for higher 
patten speeds ratios of 1.2-1.3. This is for pattern speeds from 31 to 33.5 
km s~-'^kpc~-'^, which is also out of the MW range. In this latter case, the 
effects are due to the approaching main outer resonances 4:1, 5:1 and 6:1. 
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Figure 17. Same as Fig.[T2]but for the minimum spiral density contrast. 



where substructure can be seen for the TWA. For pattern speed ra- 
tios larger than 0.8 (near corotation and high order resonances), we 
observe the same behaviour as in the PERLAS simulations. The 
comparison between both models for the minimum density con- 
trast is more drastic. The TWA model induces clear substructure 
only for ratios around 0.65-0.7 (particularly close to the 4:1 reso- 
nance). For the rest of the pattern speeds, almost no substructure is 
induced at the considered azimuths. 

Previous results referring to the pattern speeds that give im- 
portant or minor substructure do not depend substantially on the 
initial velocity dispersion or on the use of locus 1 and 2. 



10 TRANSIENT VERSUS STATIONARY EFFECTS 

A fundamental question is if spirals are long or short-lived fea- 
tures. The processes that induce spiral arms play a key role in this 
( Sellwood 2010b). With the aim to see whether the velocity dis- 
tributions can show imprints of the spiral lifetime, in this section 
we study the dependence of the structures on integration time. In 
our simulations, we have not modelled the progressive appearance 
of the spiral arms but they are included all time. In fact, some ex- 
periments showed that, given the same initial conditions, the final 
orbital structure does not change significantly with the adiabatic 
introduction of the non-axisymmetric component (Pichardo, pri- 
vate communication) or with different growth times (iMinchev et al.l 
^2010). In Fig. [TSl along columns, we show the time evolution of 
the U-V plane for 5 different representative cases (the pattern 
speed and relative spiral phase are indicated in the first panel of 
each column). We show the results for time integration intervals of 
400 Myr from to 2 Gyr, which is approximately one revolution 
period of the spiral arms. In these panels we see more separated 
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Figure 18. U-V velocity distributions 
panel of each column) for the PERLAS 



as a function of integration time for several pattern speeds or ratios VLgp/^ and different (t)sp (indicated in the first 
model with locus 1, maximum density contrast and ICl. 



groups compared to previous plots and, in general, the structures 
become separated with time. Structures in previous plots looked 
more continuous because they were composed by particles that had 
been integrated for a mix of times (between and 2 Gyr). 

The U-V plane of these simulations changes with time. This 
is due to the ongoing phase-mixing in response to spiral arm poten- 
tial. We see that the time of appearance is different for each group, 
ranging from to ^ 1200 Myr. For instance, we observe two con- 
spicuous groups in the first panel of second column (integration 
time smaller than 400 Myr). Others appear later, as the branch at 
V ^ -45 for the period 800-1200 Myr in the model of second 
column. We observe how the structures at low V or the ones at the 
extreme parts of the U-V plane require more integration time, i.e. 
a larger spiral lifetime. This is because they correspond to more 
eccentric orbits and larger radial excursions that take more time to 
reach the current region. We have also found that the time of ap- 
pearance does not depend strongly on the density contrast of the 
spiral arms. On the other hand, a hotter disc (IC2) populates some 
structures at earlier times, specially at low V . 

We see also that most structures change shape, size and posi- 
tion in the U-V plane with time. Some are formed in the central 
parts of the U-V plane and become progressively detached with 
time. Notice for example how the left structure in the third column 
shifts to negative more U. See also the branch at low V in the sec- 
ond column which becomes more populated and more extended in 
the U direction for large integration times. In general after 1200 
Myr of integration time (3-4 revolutions of the spiral arms), the 
U-V plane becomes stationary; for this and larger times all the 
groups have already appeared and preserve the same size, shape and 
position in the kinematic space. Similarly, we saw that the global 
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Figure 19. Left: Observed heliocentric U -V velocity distribution of the so- 
lar neighbourhood (data from Antoia et al. 2008). Right: Wavelet transform 
for a scale of ~ 10 km s"-*^ of the same distribution. 



density structure became stationary after integration times between 
400 and 1200 Myr (SectionO). 



11 IS IT POSSIBLE TO CONSTRAIN SPIRAL ARM 
PARAMETERS USING LOCAL STELLAR 
KINEMATICS? 

In previous sections we have explored the effects of the spiral prop- 
erties on the local velocity distribution. One would like to know 
which model parameter combination fits better the observed kine- 
matic groups in the local U-V plane. Here we discuss whether this 
constrain is currently possible. 
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Figure 20. U-V velocity distributions for several models for the PERL AS 
model with locus 1, maximum density contrast and ICl. Panels in the bot- 
tom show examples of the same model with different integration times. Or- 
ange points with error bars indicate the approximate positions of the local 
observed groups. 



11.1 Degeneracy 

The observed heliocentric velocity distribution of the solar neigh- 
bourhood is shown in Fig. [T9l ( left panel). This is obtained from the 
analysis in lAntoiaetal.l(l2008h but limited to stars with good veloc- 
ity determinations (errors in all components U, V, W ^ 2krns~^) . 
We have applied the Wavelet transform (see lAntoiaetal.ll2008h 
to this distribution with a scale of ^ lOkms"^ (right panel of 
Fig. [19]), which highlights the kinematic structures of this size. The 
main groups of this distribution are Sirius (1), Coma Berenices (2), 
Pleiades (3) and Hyades (4), groups 5 and 6 that account for the 
elongation of Hercules, and finally groups 7 and 8. To obtain ve- 
locities with respect to the LSR, we consider a solar motion of 
([/q, Vq, Wq ) = (11, 12, 7) km s-^ (Schoenrich, Binney & 
Dehnen 2010). Error bars of 5 kms~^ in this plot account for the 
uncertainty in the definition of the exact group velocity. 

We find several models, all within the established range of ev- 
idence for the MW spiral arms, that fit some of the observed groups 
in the solar neighbourhood, without any preference between these 
different "good" models. We also find that one particular group 
can be induced by different parameter combinations. We show in 
Fig. [20l several models with the PERLAS spiral arms that lead to 
kinematic groups similar to some of the observed ones. First panel 
in the first row fits Coma Berenices and presents elongations to- 
wards Hyades and groups 7 and 8. Second panel corresponds to 
the same pattern speed but a slightly different spiral phase and re- 
produces better Hyades and Sirius, showing also a small elonga- 
tion towards group 7. Third panel seems to fit Coma Berenices and 
deviations towards Pleiades. Other good fits have been found for 
the TWA model. Panels in the bottom of Fig. [20l correspond to the 
same model but different integration times. We see that especially 
for integration times from 400 to 800 Myr the model is able to 
reproduce Hyades and Sirius. This shows that the spiral lifetime 
could be also constrained. But the fact that other combinations of 
parameters have reproduced successfully these kinematic groups 
prevents us to adopt this as the best MW fitting. Notice that differ- 
ent authors have considered different observed groups and slightly 
different positions on the U-V plane due e.g. to different sam- 
ple sel ections, or d ifferent solar motion determinations such as the 
one bv lPehnen & Binnev 1998). This ambiguity, however, does not 
change the conclusions of this section and, in fact, is an additional 
limitation to constrain the model parameters. 

This degeneracy is also noticed in the literature, where dif- 
ferent models have led to the reproduction of the same group. For 
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Figure 21. Same as Fig. [12] but for a simulation including the potential 
of the Galactic bar with a pattern speed of 48 kms^-'^kpc"-'^. The outer 
Lindlbad resonance of the bar is indicated as -2:1 (B) in the right scale. 



instance Co ma Berenices have been related to the effects of the 
spiral arms ( Ouillen & M inchevll2005h and to the effects of the bar 
(Minchev et al. 2010). Our simulations have shown that the param- 
eter space is still large to obtain a unique best-fitting to the observed 
velocity distribution. 



11.2 Influence of other patterns and processes 

Apart from the mentioned current degeneracy, we should not ex- 
pect an exact equivalence between the positions of the observed 
groups and the velocity distributions from the simplified currently 
available models. We must take into account that other processes 
may have influenced the local velocity distribution: the Galactic 
bar, external processes like past accretion events, or internal disc 
processes like star formation bursts or encounters with giant molec- 
ular clouds. 

The Galactic bar is believed to infl uence considerably the lo- 
cal velocity distribution (Dehne3|2000). One may wonder how the 
velocity distributions presented in this study change under the com- 
bined action of spiral arms and bar. Several simulations have proved 
us that in such a combined case, the spiral arm imprints are still dis- 
tinguishable. Fig. [21] are velocity distributions of the same simula- 
tion as in Fig. [12] but now including also a Galactic bar. The bar is 
modelled fol l owing MW observational constraints, as explained in 
lAntoia et"aD (l2009l) . and has a pattern speed of 48kms~^kpc~^ 
and present inclination with respect to the line Sun - Galactic 
Centre of 20° . The total central mass of the m odel (bulge- f -bar) i s 
1.4 X 1O^°M0, approximately as estimated bv iPwek et°^ (Il995h . 
Of this mass, 70% belongs to the bar whereas the rest corre- 
sp onds to the bulge, with a similar proportion as the best model 
in lWeiner & Sell wood! (Il999h . This mass gives maximums for the 
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bar parameters qr and qt (defined in Section [5]) of 0.25 and 0.37, 
respectively. These values are larger than for the spiral arms but are 
achieved at inner radii. In general, we see that in the inner regions 
of the disk (up to 3 or 4 kpc) the bar dominates in strength but at 
solar radius of 8.5 kpc the strength values are comparable {qr) or 
the spiral arms are stronger {qt). 

As we see in Fig. |2T] the vast majority of structures created 
by the spiral arms are still seen in the combined case. In general, 
we see that the bar induces groups specially at the low part of 
the U-V plane. For instance, we have check that the new group 
dXV ^ — 40kms~^ in the first row appears also in the simula- 
tions with only bar. Whereas the central region of these 3 panels is 
rather similar to the spiral alone case. We find few structures that 
are not induced now (e.g. structure at low V of the third panel in 
the last row), structures that have changed their shape (third panel 
in the fourth row) or that have been shifted to a slightly differ- 
ent position of the kinematic plane (structure at low V of the first 
panel in the last row). The more severe differences are seen in the 
resona nce overlap ca se (fourth row). This is expected as, accord- 
ing to lOuillenI (l2003h . this overlap can induce widespread chaos. 
Also IChak rabartv ('2007) claimed that the potential parameter con- 
straint is not possible in these regions. But notice that we find that 
in nearby regions (^ 600 pc in galactocentric radius) the individual 
spiral effects can still be identified. To conclude, using simulations 
with spiral-only model is a valid first step to understand the iso- 
lated effects of this non-axisymmetric component on the velocity 
distribution. 

It is particularly interesting that the arms (both models 
PERLAS and TWA) can populate branch-like groups at around 

V < — 40kms~^ (Section IT2]). The V velocity of these branches 
is consistent with the V velocity of the observed Hercules struc- 
ture. This group has a V h eliocentric velocity between —40 and 
-60kms"^ (Deh nenll 19981. se e also Fig. [T9l , which with the so- 
lar motion by Sch oenrich et al. l (l201Qh . corresponds to a velocity 
with respect to the LSR of —30 to —50 km . Up to now this U- 

V plane region has been believed to be exclusively popul ated by 
the eff ects of the bar resonances. We show here and in Antoia et al.l 
(I2OO9I) that spiral and unbarred models crowd structures at these 
negative V. Although its shape is not exactly equal to the Hercules 
branch, especially regarding its observed inclination in the U-V 
plane and average radial motion U < 0, we have just shown that the 
Galactic bar action can modify and shift the spiral-induced groups 
to more negative [/. As a ba r-only model is also able to induce a 
grou p similar to Hercules (Dehnenl bood : IFuxI|2001I : lAntoia et al.' 
I2OO9) . it is difficult to favour the bar or the combination of spiral 
arms and bar as the cause of this group. 

The characteristics of the Galactic bar are still quite uncertain. 
For instance, determinations of the bar's an gle with respect to the 
line Sun - Galactic Centre range from 14° ( Freudenreich "1998*) to 
40 - 45° (iHammerslev et al.ll2000l : iBeniamin et al. 2005). More- 
over, the existence of one or two bars in the MW is currently being 
debated. It would be worth exploring in detail the kinematic ef- 
fects of the combined spiral-bar allowed parameter space, although 
this is out of the scope of the present study. However, note that a 
different bar's orientation could produce some different kinematic 
groups, but will not change the conclusions of the section. This is 
not strictly true if we change the relative strength of the bar with 
respect to the arms or their pattern speeds. These may lead to kine- 
matic distributions dominated by one of the patterns. Two extremes 
cases would be when one pattern exceeds the other one in strength 
or when its resonances are much closer to the considered regions. 
However, the ranges for the parameters of the spiral arms and bar 



in our models based on literature seem to indicate that the MW is 
far from these regimes. Simulations with other bar parameters, not 
presented in this study, have confirmed that the results of this sec- 
tion do not depend strongly on the parameters provided they are in 
the believed MW range. 

We could also consider including other kind of patterns in the 
disc as proposed in the literature. For instance, several models have 
shown that galaxies might have more than one spiral arm mode 
(with differ ent patter n speeds) at different radius or overlapping in 
radius (e.g. lRautiainen & Salo 1999), or with pattern speed chang- 
ing with radius (e.g. Dobbs..201Q) . Some of the modes could be 
corotating with the bar at the bar's end, or coupled with the bar 
through their respective resonances. These possibilities are being 
explored for other galaxies (Meidt, Rand & Merrifield 2009) but 
have hardly been addressed for the MW. These must be future 
points to take into account in our modelling. 



12 SUMMARY AND CONCLUSIONS 

Few studies have focused on the spiral arm effects on the local ve- 
locity distribution. Moreover, all of them modelled the spiral arm 
potential following the TWA. We have seen that, despite being un- 
certain, the observational evidences for the MW arms density con- 
trast, and especially, the pitch angle, suggest that the assumptions 
for self-consistency of the TWA in the MW case are, at least, doubt- 
fully satisfied. Here we have studied the spiral arm effects on the 
kinematics of the solar neighbourhood with the TWA but also with 
a different spiral arm model, the PERLAS model. While in the 
TWA, the spiral arms are a small perturbative term of the poten- 
tial, modelled as a cosine function, the PERLAS spiral arms cor- 
respond to an independent 3D mass distribution, from which the 
gravitational forces are derived. We have seen that both radial and 
tangential forces of this model present more abrupt features and 
different disc locations for the minima, maxima and points. 

Here we aimed to determine the conditions that favour the ap- 
pearance of kinematic groups, such as the ones that we observe in 
the solar neighbourhood. To do this, first we have tuned these two 
models to the latest observational determinations for the properties 
of the MW spiral arms. For many of these properties only a certain 
range of evidence could be established. Next, we have performed 
test particle simulations with both models considering these ranges, 
different initial conditions and integration times. Our analysis indi- 
cates that: 

(i) A significant part of the allowed parameter space, especially 
for the PERLAS model, favours the triggering of kinematic groups 
with different shapes, sizes and inclinations, such as the observed 
ones in the solar neighbourhood. This shows that it is feasible that 
some of the observed moving groups have a dynamical origin re- 
lated with the spiral arms. 

(ii) The velocity distribution is certainly sensitive to the relative 
spiral arm phase and, especially, to the pattern speed. Changes in 
the kinematic groups are significant if one moves only ^ 0.6 kpc in 
radius. But the U-V plane changes more slowly with azimuth and 
^ 2 kpc are needed in this direction to detect important differences. 
Due to this, the pattern speed could be better constrained using the 
observed kinematic groups (within an error of < 2kms"^kpc"^) 
than the relative spiral phase (with a precision of ~ 15°). However, 
one would need to break the degeneracy mentioned below by other 
means such as reducing the free parameter space. 

(iii) For both models and for all density contrasts, within the ob- 
servational MW range, the spiral arms induce strong imprints for 
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pattern speeds around 17kms~^kpc~^. This corresponds to re- 
gions close to the 4:1 inner resonance. No substructure at all is in- 
duced close to corotation or high order resonances (m > 6), which 
corresponds to pattern speeds of 20.5 to 30 kms~^kpc~^. 

(iv) Changes in spiral strength produce no significant differ- 
ences in most cases, which makes difficult to constrain this param- 
eter. Some groups in the extremes of the U-V plane are only pop- 
ulated for the maximum density contrast case, which could help to 
establish the spiral strength. The effects of a different pitch angle 
seem difficult to differentiate given the observational range for this 
parameter. 

(v) The spiral arms induce groups in a large region of the U-V 
plane, including the low V part. For instance, the arms, and not only 
the bar, can crowd the region of the Hercules group. Previous stud- 
ies associated the spiral arm influence mainly to the central parts of 
the U-V plane, and the bar to low angular momentum structures, 
such as Hercules (Guillen & Minchev 2005; Dehnen 2000). 

(vi) The kinematic groups that are dynamically induced by the 
spiral arms depend on integration time. The time of appearance is 
different for each group, ranging from to 1200 Myr. The struc- 
tures at the extreme parts of the U-V plane require more integra- 
tion time, i.e. a larger spiral lifetime. The early appearance of some 
groups demonstrates that even recent spiral arms (< 400 Myr) 
may produce strong kinematic structures. Structures also change 
shape, size and position in the U-V plane until a maximum time of 
^ 1200 Myr, when the U-V plane becomes stationary. 

(vii) Each of the structures seen in the U-V plane is composed 
by particles with a wide range of integration times and a different 
minimum time. This is encouraging as the study of the stellar ages 
in the observed moving groups also reveal similar characteristics. 
For instance, see evolution in the Sirius and Hyades groups partic- 
ularly for stars younger than 2 Gyr (fig. 13 in lAntoia et al.l2 008) or 
the wide age distri bution and differ ent minimum age in a particular 
group (fig. 14 in Antoia et al. ll2008h . 

(viii) Models obtained from our simulations with the spiral arms 
that reproduce the local observed velocity distribution are degener- 
ated. Groups such as the ones that are observed in the solar neigh- 
bourhood can be induced by different model parameter combina- 
tions. For instance, several models create groups such as Hyades, 
Pleiades or Sirius. 

(ix) In most of our simulations where both the spiral arms and 
the bar are included, individual imprints of the bar and the arms 
can still be identified in the final velocity distributions. This means 
that using simulations with spiral-only model is a valid first step to 
understand the isolated kinematic effects of this non-axisymmetric 
component. 

(x) The stellar response near solar positions of the TWA spi- 
ral arms and the PERLAS model is significantly different, both in 
configuration space and kinematics. Both models are able to in- 
duce several kinematic groups but the velocity distributions and the 
groups are rather different. 

(xi) The PERLAS model gives more substructure than the TWA 
given the same density contrast, and more substructure for a larger 
pattern speed range or, equivalently, a larger radius range. For the 
PERLAS model and both the maximum and minimum density 
contrast, we see clear and rich kinematic substructure near solar 
azimuth for pattern speeds from 15 to 19.5 kms~^kpc~^ (regions 
near the inner 3: 1, 4: 1 and 5:1 resonances). By contrast, this pattern 
speed range is smaller for the TWA model even in the high density 
contrast case. On the other hand, for the minimum density contrast, 
the TWA model induces some substructure only for a narrow range 
of pattern speeds (around 16.7-18 km s~^kpc~^). 



To conclude, the spiral arms induce strong imprints on the ve- 
locity distributions at solar regions and that these are sensitive to 
some arm properties. These indicate that kinematics can be used as 
one of the constraints on the current uncertainty about spiral struc- 
ture of our Galaxy, e.g., in the pattern speed, strength, orientation 
and lifetime. 

The local distribution function is significantly sensitive to the 
used model, even if they are adjusted to reproduce the same obser- 
vational constraints for the spiral structure of the MW. This stresses 
the importance of the specific spiral arm gravitational potential 
modelling for the MW studies, as this can induce significant devi- 
ations on their orbital dynamics. It is promising that the kinematic 
structures can help to improve this modelling or to test whether the 
MW spiral arms are weak and tightly wound, following the TWA. 

The constraint of the MW spiral structure properties, how- 
ever, is presently not straightforward. This is due, first, to the men- 
tioned degeneracy. Our simulations have shown that the parame- 
ters space is still too large to obtain a unique fitting to the ob- 
served velocity distribution. Second, the uncertainty in the solar 
motion propagates to the velocities of the observed groups, mak- 
ing the fitting between observations and simulations imprecise. The 
problem is complex as, precisely, the presence of kinematic sub- 
structure complicates the solar motion determination. Third, for 
this constraint we should take into account the other processes that 
have influenced and sculpted the real velocity distributions (other 
non-axisymmetries, external accretion effects, dispersion by giant 
molecular clouds,etc). As an example, we have identified effects 
due to the combination of bar and spiral arms. For instance, struc- 
tures that appear centred in the U axis in the spiral-only models are 
shifted to negative U due to the combined action, inducing a group 
similar to the observed Hercules structure. Nevertheless, none of 
the existing current Galaxy models is complex enough to include 
all these processes at the same time. 

To break the degeneracy and use the kinematic groups to con- 
strain the MW large scale structure, a smaller parameter space for 
the spiral properties and/or data from velocity distributions at dif- 
ferent positions of the MW disc are needed. With the increase of 
knowledge on the MW structure and evolution, we will have to 
cope with a more complex scenario and a variety of processes that 
can play a role in the formation of moving groups. Future models 
will be more realistic but more complex to interpret. Besides, while 
now the spatial study of the observed moving groups is limited by 
the extension and precision of the current observational samples, 
this will soon change with the advent of data from new surveys 
(USNO, UCAC3, 2MASS, RAVE, SEGUE, PanStars, LAMOST, 
Gaia). Explorations of the spiral arm effects like the present study 
will help us to interpret more sophisticated future models as well as 
velocity distributions in nearby regions coming from surveys such 
as Gaia. 



ACKNOWLEDGMENTS 

We thank the referee, P. Rautiainen, for useful suggestions that 
helped improve this manuscript. We acknowledge funding sup- 
port from the European Research Council under ERC-StG grant 
GALACTICA-24027, MICINN (Spanish Ministry of Science and 
Innovation) - FEDER through grant AYA2009-14648-C02-01 and 
CONSOLIDER CSD2007-00050. Some of the simulations were 
run at the HP CP 4000 cluster (KanBalam) in the DGSCA/UNAM. 



18 Antojaetal. 



REFERENCES 

Allen C, Santillan A., 1991, Rev. Mexicana Astron. Astrofis., 22, 
255 

Antoja T., Figueras R, Fernandez, D., Torra J., 2008, A&A, 490, 
135 

Antoja T., Valenzuela O., Pichardo B., Moreno E., Figueras D., 

Fernandez, D., 2009, ApJ, 700, 78 
Antoja T., Figueras R, Torra J., Valenzuela O., Pichardo B., 2010, 

Lecture Notes and Essays in Astrophys. IV, 4, 13 
Athanassoula E., Bienayme O., Martinet & L., Pfenniger D., 

1983, A&A, 127, 349 
Aumer M., Binney JJ., 2009, MNRAS, 397, 1286 
Benjamin R.A. et al., 2005, ApJ, 630, 149 
Bensby T., Oey M.S., Feltzing S., Gustafsson B., 2007, ApJ, 655, 

89 

Binney J., Tremaine S., 2008, Galactic Dynamics, Princeton Uni- 
versity Press Princeton, New Jersey 

Block, D.L., Witt, A.N., Grosb0l,P., Stockton, A., Moneti, A. 
1994, A&A, 288, 383 

Block D.L,, Buta R., Knapen J.H, Elmegreen D.M., Elmegreen B. 
G., Puerari L, 2004, AJ, 128, 183 

Chakrabarty D., 2007, A&A, 467, 145 

Churchwell E. et al, 2009, PASP, 121, 213 

Combes, R, Sanders, R. H. 1981, A&A, 96, 164 

Contopoulos G., Grosb0l R, 1986, A&A, 155, 11 

Dehnen W., 1998, AJ, 115, 2384 

Dehnen W., 1999, AJ, 118, 1201 

Dehnen W., 2000, AJ, 119, 800 

Dehnen W., Binney J., 1998, MNRAS, 298, 387 

De Simone R.S., Wu X., Tremaine S., 2004, MNRAS, 350, 627 

Drimmel R., Spergel D.N., 2001, ApJ, 556, 181 

Dobbs, C.L., 2010, preprint (astro-ph/1002.1269D) 

Dwek, E., et al., 1995, ApJ, 445, 716 

Elmegreen D.M., Elmegreen B.G., 1990, ApJ, 355, 52 

Eskridge, RB., et al., 2002, ApJS, 143, 73E 

Franco J., Martos M., Pichardo B., Jongsoo K., 2002, ASPC, 275, 
343 

Freudenreich H.T., 1998 ApJ, 492, 495 
Fux R., 2001, ApJ, 373,511 

Georgelin Y.M., Georgelin Y.P, 1976, A&A, 49, 57 

Gerhard O., 2010, preprint (astro-ph/1003.2489G) 

Gerola H., Seiden RE., 1978, ApJ, 223, 129 

Grosb0l R, Patsis RA., Pompei E., 2004, A&AS, 423, 849 

Hammersley, P.L., Garzon, R, Mahoney, T. J., Lopez-Corredoira, 

M., Torres, M.A.R, 2000, MNRAS, 317, 45 
Hemquist L., 1993, ApJS, 86, 389 

Holmberg J., Nordstrom B., Andersen J., 2009, A&A, 501, 941 
Kalnajs A. J., 1991, Dynamics of Disk Galaxies, eds. B. Sun- 
delius, 323 

Kendall S., Kennicutt R. C, Clarke C, 2011, preprint (astro- 

ph/1 101.5764) 
Kennicutt R.C., Hodge R, 1982, ApJ, 253, 101 
Kent S.M., 1992, ApJ, 387, 181 
Kim W.-T, 2009, lAUS, 254, 313 
Lin C.C., Yuan C, Shu RH., 1969, ApJ, 155, 721 
Martos M.A., Cox D.R, 1998, ApJ, 509, 703 
Martos M., Yanez M., Hernandez X., Moreno E., Pichardo B., 

2004, JKAS, 37, 199 
Meidt, S.E., Rand, R.J., Merrifield, M.R., 2009, ApJ, 702, 277 
Minchev L, Boily C, Siebert A., Bienayme O., 2010, MNRAS, 

407, 2122 



Miyamoto M., Nagai R., PASJ, 27, 533 
Oort J.H., Kerr RJ., Westerhout G., 1958, MNRAS, 118, 379 
Patsis RA., Contopoulos G., Grosb0l R, 1991, A&A, 243, 373 
Pichardo B., Martos M., Moreno E., Espresate J., 2003, ApJ, 582, 
230 

Press W.H., Teukolsky S.A., Vetterling W.T., Flannery B.R, 1992, 
Numeric al Recipes in Fortran 77: The Art of Scientific Comput- 
ing, 2nd. ed. (Cambridge University Press: Cambridge) 

Quillen A.C., 2003, AJ, 125, 785 

Quillen A.C., & Minchev L, 2005, ApJ, 130,576 

Rautiainen, R, Salo, H., 1999, A&A, 348, 737 

Reid M.J. et al., 2009, ApJ, 700, 137 

Rix H.-W., Zaritsky D., 1995, ApJ, 447, 82 

Roberts W.W., Jr., Huntley J.M., van Albada G.D., 1979, ApJ, 
233, 67 

Romero-Gomez M., Athanassoula E., Masdemont J.J., Garcia- 

Gomez C, 2007, A&AS, 472, 6 
Schoenrich R., Binney J., Dehnen W., 2010, MNRAS, 403, 1829 
Schmidt M., 1956, Bull. Astron. Inst. Neth. 13, 468, 15 
Sellwood J.A., 2010a, MNRAS, 409, 145 
Sellwood J.A., 2010b, MNRAS, 410, 1637 
Simonson S.C., 1970, A&AS, 9, 163 

Skuljan J., Heamshaw J. B., Cottrell RL., 1999, MNRAS, 308, 
731 

Taylor J.H., Cordes J.M., 1993, ApJ, 411, 674 

Toomre A., 1981, in The structure and evolution of normal galax- 
ies; Proceedings of the Advanced Study Institute, Cambridge 
University Press, Cambridge , p. 1 1 1 

Vallee J.R, 2008, ApJ, 135, 1301 

VogHs N., Tsoutsis R, Efthymiopoulos C, 2006, MNRAS, 373, 
280 

Weiner, B.J., Sellwood, J.A., 1999, ApJ, 524, 112 



